Ensemble inequivalence, bicritical points and azeotropy for generalized Fofonoff flows 
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We present a theoretical description for the equilibrium states of a large class of models of two- 
dimensional and geophysical flows, in arbitrary domains. We account for the existence of ensemble 
inequivalence and negative specific heat in those models, for the first time using explicit computa- 
tions. We give exact theoretical computation of a criteria to determine phase transition location 
and type. Strikingly, this criteria does not depend on the model, but only on the domain geometry. 
We report the first example of bicritical points and second order azeotropy in the context of systems 
with long range interactions. 
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In many fields of physics, the particle or fields dynam- 
ics is not governed by local interactions. For instance for 
self gravitating stars in astrophysics [l|, for vortices 
in two dimensional and geophysical flows [3, 0, [3], for 
unscreened plasma or models describing interactions be- 
tween waves and particles @ , the interaction potential is 
not integrable 0|. Recently, a new light was shed on the 
equilibrium statistical mechanics of such systems with 
long range interactions : there has been a mathematical 
characterization of ensemble inequivalence @ , a study of 
several simple models [3, Q3] , and a full classification of 
phase transitions and of ensemble inequivalence [ll| . 

One of the promising field of application for the statis- 
tical mechanics of systems with long range interactions, is 
the statistical prediction of large scale geophysical flows. 
For instance, the structure of Jupiter's troposphere has 
been successfully explained using the Robert-Sommeria- 
Miller (RSM) equilibrium theory [H] [13]. One of the 
major scope of this field is to go towards earth ocean 
applications. All textbook in oceanography present the 
Fofonoff flows which have played an important historical 
role in that field [3]. In this letter, we propose a theoret- 
ical description of such flows in the context of the statisti- 
cal theories which, for the first time, relates its properties 
to phase transitions (see Fig. [TJ), negative specific heat 
and ensemble inequivalence. 

One of the striking features of the equilibrium theory 
of systems with long range interactions is the generic ex- 
istence of negative specific heat. This strange phenomena 
is possible as a consequence of the lack of additivity of 
the energy and is related to the inequivalence between 
the microcanonical and canonical ensemble of statisti- 
cal physics. This was first predicted in the context of 
astrophysics [13]. For two dimensional flows, existence 
of such inequivalence has been matematically proven for 
point vortices [l6| (without explicit computation), and 
numerically observed in a particular situation of a Quasi- 
Geostrophic (QG) model [13], and in a Monte Carlo 
study of points vortices in a disk [lij [. One of the novelty 
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Figure 1: Geometry governed phase transition : second order 
phase transition from Fofonoff modes in a domain with axial 
symmetry to dipole solutions, breaking the symmetry, only 
when the domain is sufficiently stretched (r > re) ■ 



and achievement of the current work, is that we predict 
such ensemble inequivalence, with exact theoretical com- 
putation of the associated phase transitions, for a very 
large class of models including the Euler equation or QG 
models. 

In the context of systems with long range interactions, 
a classification of phase transitions associated to ensem- 
ble inequivalence has been proposed 
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transitions predicted have never been observed, neither 
in models, nor in real physical systems. One of the main 
interest of the current work, is the finding of two ex- 
amples of such unobserved phase transitions : bicritical 
points (a bifurcation from a first order phase transition 
towards two second order phase transitions) and second 
order azeotropy (the simultaneous appearance of two sec- 
ond order phase transitions at a bifurcation). We prove 
that those phase transitions are governed by the geom- 
etry of the domain in which the flow takes place. We 
explain how to easily compute the transition point for 
any domain geometry. 

Euler and QG equations This letter describes phase 
transitions existing in a broad ensemble of models : 2D 
Euler flows, one layer QG models in a closed domain 
T>. The common character of all of those models comes 
from the fact that they can all be expressed as a quasi- 
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2D transport equation dtq + u • Vg = 0. For the one 
layer QG model, the potential vorticity (PV) q = Atp — 
tp/R 2 + h is a scalar ; h(x,y) is the topography and R 
is the Rossby radius of deformation. The velocity field 
is related to q via tp : u = e 2 x 'Vip. The case h = 0, 
R = +oo corresponds to the Euler equation. 

For all of these models we use an impermeability 
condition at the boundaries. The energy of such sys- 
tems can always be written £[q] = — \{0[q — h],q— h) 
where O is a symmetric linear operator. For instance, 
£ = \{{Vi>f + tJj 2 /R 2 ) in the 1-1/2 layer QG model, 
which corresponds to Q = (A — i?~ 2 ) _1 . In all these 
models, both energy and circulation C[q] = {q) are con- 
served quantities. 

Variational problem In the following we consider the 
solutions of the variational problem : 



S(E,T) = max{S[ 9 ] 



£[q] =Ek C[q] = T} , (1) 



where S is the entropy of the PV field q : S[q] = (s(q)), 
with s a concave function, and S (E, T) is the equilibrium 
entropy. Such a variational problem may be interpreted 
in two ways. First, in the Robert-Sommeria-Miller sta- 
tistical mechanics of the Euler equation, one obtain a 
much more complex variational problem which involves 
the usual Maxwell-Boltzmann entropy and which is con- 
strained on the whole initial vorticity distribution. One 
can prove that if this vorticity distribution is treated 
canonically, then the statistical equilibrium verify the 
variational problem (TT1) 1 191 . Moreover any solution to ([1]) 
is a RSM equilibrium [13]. An alternative interpretation 
is to assume directly that, in some physical situations, 
one has to consider the most probable state with respect 
to a prior vorticity distribution, which can be related to 
the shape of s [§]. 

To compute critical points of the variational problem 
([l]), we introduce two Lagrange parameters (3 and 7 as- 
sociated respectively with the energy and the circulation 
conservation. These critical points are stationary solu- 
tions for the initial transport equation : q = f(ip), with 
/ (ip) = s'^ 1 (—/3ip + 7). In all of the following, we study 
the case of a quadratic entropy S[q] = ( — g 2 /2). This 
leads to a linear vorticity-stream function relationship 
(/ ("0) = Pi 1 — l)- The original Fofonoff solution corre- 
sponds to the particular case (3 3> 1 and h(x,y) = y. 
Among the states we study, only the ones with /? lower 
than the phase transition value have already been de- 
scribed as Fofonoff flows or in the context of the Kraich- 
nan statistical mechanics 2(], 21 1 . In the context of the 
RSM statistical mechanics, in the case of Euler equation, 
Chavanis and Sommeria [22| found a criteria for the ex- 
istence of a transition from a monopole to a dipole when 
increasing the energy (as in figure]!]). By using a different 
method (by solving directly the variational problem) , we 
generalize these results to a wide class of model, analyze 
for the first time ensemble inequivalence, find unobserved 



phase transitions and establish the relations of these to 
Fofonoff flows. 

Dual variational problems Our problem is to find the 
minimum of a quadratic functional <JTJ> , taking into ac- 
count the constraints on circulation (linear) and energy 
(quadratic with possibly a linear contribution). This 
will be referred as the microcanonical problem. Deal- 
ing with unconstrained variational problems is much eas- 
ier than dealing with constrained ones: solutions for a 
variational problem are necessarily solutions for a more 
constrained dual problem @|. Moreover, when all the 
possible constraint values are achieved in the less con- 
straint ensemble, we are in a situation of ensemble equiv- 
alence. The study of the unconstrained variational prob- 



lems is then sufficient [ll|. We will thus consider, by 
relaxing one or both constraints, two dual problems : a) 
canonical, by relaxing the energy constraint, with the 
free energy F{fi, T) = min 9 {-%] + £[q] \ C[q] = T}; 
b) grand canonical, by relaxing both energy and circu- 
lation constraints, with the thermodynamical potential 
J(/3, 7) = min 9 {— S[q] + f3 £[q] + jC [q] } . It is natural 
to consider first the grand canonical ensemble. If en- 
semble inequivalence does exist, we will then study a 
more constrained variational problem, and so on, until 
the whole range of E and T has been covered. Notice 
that in our case, the ensemble of accessible values for 
(£, C) is the half plane E > 0. 

Solutions of quadratic variational problems For all 
variational problem to be considered in the following, we 
will look for the minimum of a quadratic functional, with 
a linear part. Let us call Q the purely quadratic part and 
L the linear part of this functionnal. Then we have three 
cases 

1. The smallest eigenvalue of Q is strictly positive. The 
minimum exists and is achieved by a unique minimizer. 

2. At least one eigenvalue of Q is strictly negative. There 
is no minimum (the infimum is —00). 

3. The smallest eigenvalue of Q is zero (with eigenfunction 
eo). If Leo = (case 3a), the maximum exists and each 
state of the neutral direction {aeo} is a minimizer. If 
Leo 7^ 0, (case 3b) then no minimum exist. 

The grand canonical ensemble We consider in this 
part h = 0, so the energy is purely quadratic. We look 
for the minimum of = — S + (3 £ + 7C We intro- 
duce the projections qi of the PV q on a complete or- 
thonormal basis of eigenfunctions ei(x, y) of the operator 
O (0[ei] = ej/A,) (see the definition of the energy) and 
find: 

j[q\ = J2 s v i 1 - pi x *) w + J2 

i.j>l i>l 

where Sij is the Kronecker symbol, and where the Aj 
(all negative) are in decreasing order. One can see that 
the quadratic part is diagonal, and that all eigenvalues 
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are strictly positive if and only if (3 > Ai(case 1. above) 
. If (3 — Ai(case 3), a neutral direction exists if and 
only if 7 = (case 3a). Thus, grand canonical solu- 
tions exist only for (3 > Ai or (/3 = Ai and 7 = 0). By 
computing the energy and circulation of all those states, 
we prove that it exists a unique solution at each point 
in the diagram (E,T), below a parabola V of equation 

E = — T 2 Ai/ ^2(ei) 2 ^. Because the values of energies 
above the parabola V are not reached, we conclude that 
there is ensemble inequivalence for parameters in this 
area. As explained previously, we have to consider a more 
constrained variational problem to find solutions in this 
area. 

The canonical ensemble We now look for the mini- 
mum of —S + (3 £ for a fixed circulation V. In order to 
take into account this constraint, we express one coordi- 
nate in term of the others : q\ = (T — J2i1i( e i)) I ( e i)- 
Then, we have to minimize 



i,j>2 



i>2 



(ei> 2 V Aj' 



A J < ei > 2 )<Mj 



without constraint. 

The linear operator Q, associated to the quadratic part 
of J 7 , is not diagonal in the basis {e^}. We first notice that 
if the domain geometry admits one or more symmetries, 
it generically exists eigenfunctions having the property 
(ei) = 0. In the subspace spanned by all those eigenfunc- 
tions, Q is diagonal, and its smallest eigenvalue is posi- 
tive as long as (3 > P®, where /3j is the greatest Ai on this 
subspace. Then we look for the value of /? such that the 
smallest eigenvalue of Q is zero in the subspace spanned 
by eigenfunctions with (ej) ^ 0. Let us call (3* this value, 
and q* the corresponding eigenfunction : Q[q* c ] = 0. Us- 
ing this last equation, we prove that (3* is the greatest 
zero of the function f(x) = l — xJ2i>i( e i) 2 /( x ~^i) ■ We 
conclude that there is a single solution to the variational 
problem for (3 > max(/3j ) , /?*) (case 1) and no solution 
for (3 < max (/3?, (3* c ) (case 2). When (3 = max /?*), 
to discuss the existence of a neutral direction, we distin- 
guish two cases according to the sign of /3° — /?* : 

i) 0( < /3* we then consider /? = /?*. We are in case 3a for 

T = and in case 3b for Y / 0. 

ii) /3i > P* we then consider (3 = (3°, We are in case 3a 

whatever the value of T. 

In case i), if T ^ 0, all energy value are reached in the 
canonical ensemble. If T = 0, the solutions of the neutral 
directions are aq*. Thus, varying a, all energy values are 
reached ; there is two canonical solutions corresponding 
to each energy E, depending on the sign of a. In case ii), 



YiYHfj^po E(0) — E-p (T) oc r 2 . This defines a parabola 
Vq. Then, whatever the value of T, there is a unique 
canonical solution q (/3, T) for each point of the diagram 
(E,T) below Vq. The canonical solutions of the neutral 
directions is q = ae\ + lim^^o q ((3, T) where e\ is the 
eigenfunction associated with 13®. Varying a, all energy 
larger than E-p (T) are reached. For each energy above 
the parabola Vq, it exists two canonical solutions, de- 
pending on the sign of a. In both cases we find that 
all circulation and energy values have been reached by 
canonical solutions. We conclude that microcanonical 
and canonical ensembles are equivalent. 

Description of phase transitions. In case i), at 
fixed energy, we can show by a direct computation that 
7(T) = ff is discontinuous in T = 0. This means that 
there is a microcanonical first order transition (see point 
G figure (21-c) . In case ii), one can show that §^ is dis- 
continuous when (E 7 T) belongs to Vo (see points B and 
D, figure [2]-b) . The ensemble inequivalence area is as- 
sociated with the existence of a first order transition in 
the ensemble with only one constraint on the energy (see 
the corresponding Maxwell constructions on figure[2]-b,c). 
Similarly, if we now fix the circulation, there is a discon- 
tinuity of when (E, T) belongs to Vq. It means that 
■Po is a line of second order phase transition for canonical 
and microcanonical ensembles. 

General criteria. The main interest of the abstract 
previous analysis is to conclude that all of the models 
considered will behave according to only two types of 
phase diagram structures. The difference between the 
two classes of systems is the existence of either first or sec- 
ond order phase transitions, corresponding respectively 
to case i) and ii) . If the domain has no symmetry, only 
case i) is possible (generically). If there is a symme- 
try, both cases are possible. The criteria for class i) or 
class ii) systems is the sign of (3® — (3*, that can be eas- 
ily computed for a system at hand. The same criteria 
was obtained in [22j for the Euler equation, using a dif- 
ferent method. Very interestingly, this criteria does not 
depend on the model considered (it does not depend on 
the Rossby radius of deformation R), but only on the 
domain shape. For some classes of domain geometries, 
as ellipses for instance, we are always in case ii). More 
generally, any domain geometry sufficiently stretched in a 
direction perpendicular to its symmetry axis is in case ii). 
This is for instance the case of a rectangular domain The 
transition from systems of type i) to systems of type ii), 
when the geometry is modified, leads to very interesting 
phenomena that are described now. 

The bicritical point On figure O we consider a fixed 
energy, and present the phase diagram in the (T, t) plane, 
where r is a parameter characterizing the aspect ratio 
of the domain (horizontal over vertical width). In the 
microcanonical ensemble, there is a bifurcation from a 
first order transition line to two second order transition 
lines at a critical value r = t c . Such a bifurcation is 
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Figure 2: First observation of a bicritical point (change from a 
first order to two second order phase transitions) in a system 
with long range interactions. Euler equation in a rectangular 
domain of aspect ratio r. Green line : second order phase 
transition (discontinuity of dj/dF in b). Red line : 1st order 
phase transition (discontinuity of 7 in c). Insets are projection 
of the Entropy S [q] in a plane (5? , qi ) for fixed energy and 
circulation. 
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Figure 3: First observation of second order azeotropy (at en- 
ergy Ea there is simultaneous appearance of two second order 
phase transitions), a) dj/dV is continuous ; b and c) d'y/dY 
is not continuous in Fi and F2. 



referred as a bicritical point (see [lll|). 

With a topography We now give the main striking fea- 
tures arising when adding the term h, especially the new 
unobserved phase transitions. Concerning the existence 
of first and second order phase transitions in the micro- 
canonical ensemble, there is now three possibilities. If 
— P* < 0, the phase diagram is similar to the one of 
case i). If — /3* > and (h, e°) = 0, then we are in the 
case ii) , except that the minima of V and Vo are no more 
at the same place in the diagram (E,T), If 0% — 0* > 
and (h, ej) ^ 0, there is neither second nor first order 
transition. 

Second order Azeotropy In case ii) with topography, 
if we consider the energy as an external parameter, there 
is the simultaneous appearance of two second order phase 
transitions in the microcanonical ensemble (see figure [3] 
and the corresponding flows figure Q]), which is the signa- 
ture of second order azeotropy. 

Conclusion We report the generic existence of ensem- 
ble inequivalence and of new phase transitions in a large 
class of 2D flows. All phase transitions presented here 
appear in the inequivalence ensemble area. They are in 
that respect a signature of such an inequivalence in a 



long range interacting system, which has never been ob- 
served experimentally. The observation of those transi- 
tions could be carried in laboratory experiments on quasi 
2D flows. This could be done by using either magnetized 
electron columns [13] or three dimensionnal tanks with 
small height compared to the horizontal scale, with a 
further ordering (strong rotation or a transverse mag- 
netic field). The interest for ocean applications will be 
described in a companion paper, as well as the detailed 
complete computations, generalization to multi level QG 
equations and analysis of the ensemble inequivalence as- 
sociated to the bicritical points. 
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